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Abstract 

High resolution and seeing-free spectroscopic observation of a decaying sunspot was done with the Solar 
Optical Telescope aboard Hinode satellite. The target was NOAA 10944 located in the west side of the 
solar surface from March 2 to March 4, 2007. The umbra included many umbral dots (UDs) with size of 
—300 km in continuum light. We report the magnetic structures and Doppler velocity fields around UDs, 
based on the Milne-Eddington inversion of the two iron absorption lines at 6302A. 

The histograms of magnetic field strength(5), inclination angle (i), and Doppler velocity (u) of UDs 
showed a center-to-limb variation. Observed at disk center, UDs had (i)slightly smaller field strength 
(AB = —17 Gauss) and (ii)rclative blue shifts (Av =28 m s _1 ) compared to their surroundings. When the 
sunspot got close to the limb, UDs and their surroundings showed almost no difference in the magnetic 
and Doppler values. This center-to- limb variation can be understood by the formation height difference 
in a cusp-shaped magnetized atmosphere around UDs, due to the weakly magnetized hot gas intrusion. 
In addition, some UDs showed oscillatory light curves with multiple peaks around 10 min, which may 
indicate the presence of the oscillatory convection. We discuss our results in the frameworks of two 
theoretical models, the monolithic model (Schussler & Vogler 2006) and the field-free intrusion model 
(Spruit & Scharmer 2006). 
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1. Introduction 

The sunspot is one of the most prominent structures 
in the solar photosphere, although there are many re- 
lated unsolved problems remaining even today. One of 
them is the source of energy transport in sunspots. It 
is known that the radiative energy alone is not sufficient 
for accounting for the observed brightness of sunspots, so 
another form of convectivc energy transport is necessary 
(Deinzcr 1965). The study of umbral dots (UDs), which 
are tiny bright points in the umbra, are essential for un- 
derstanding the energy transport in sunspots, since UDs 
are considered to be a manifestation of the convection. 

Parker (1979) suggested in his "spaghetti" model, that 
UDs are the radiative signatures of the top parts of field- 
free convectivc plumes. The field-free plumes intrude from 
below the visible surface into a gap between magnetic bun- 
dles in a cluster- type sunspot. These plumes are accom- 
panied by smaller magnetic field strength, substantial up- 
flows within UDs, and a cusp-shaped magnetic structure 
(Spruit & Scharmer 2006). Another promising mechanism 
is the magneto-convection in a monolithic sunspot. The 
monolithic model considers a sunspot as the aggregation 
of uniform vertically thin columns, and UDs as a natural 
result of the overstable oscillatory convection, which is a 
preferred mode just below the photosphere (Weiss et al. 
2002; Schussler & Vogler 2006). The monolithic model 
predicts smaller field strength, upflows in the center of 



UDs in addition to downflows at their boundaries, and a 
cusp-like structure. 

Few spectroscopic observations of UDs have been done 
so far, because of their tiny size (less than 0"5) and low 
brightness. Only in recent days, some spectroscopic works 
have been published (Pahlke & Wiehr 1990; Lites et al. 
1991). Wiehr & Degenhardt (1993) observed UDs in the 
lines Fe 6843A and Ca 6103A, and found field strength re- 
duction up to 20% and flatter field inclination only in the 
lower layer. Weaker field (—500 Gauss) with more hori- 
zontal orientation (-10°) in UDs was reported in Socas- 
Navarro et al. (2004). 

As for the Doppler velocity field, substantial upward ve- 
locity (— 1 km s _1 ) is observed in the lower photosphere. 
However in the lines formed in the upper photosphere, no 
strong velocity field is related to individual UDs. Rimmclc 
(2004) found upflows in excess of 1 km s" 1 in C I 5380A 
(lower photosphere) line, while no strong upflow in Fc 
I 5576A (upper photosphere) line. Socas-Navarro et al. 
(2004) also found upflows of -250 m s" 1 in Fe 6303.46A 
(lower photosphere) line, while no obvious upflow in Fe 
6302. 5A (upper photosphere) line. Recently, Bharti et al. 
(2007a) found upward velocity on the order of 400 m s -1 , 
surrounded by narrow downflow regions with —300 m s _1 
in Fe I 5 5 76 A line. 

The spectro-polarimeter (SP) on board of the Hinode 
Solar Optical Telescope (SOT) (Tsuneta et al. 2008; 
Suematsu et al. 2008; Ichimoto et al. 2008; Shimizu et al. 
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2008) made it possible to observe diffuse UDs at the center 
of the umbra, with resolution limit of 0"3 in a highly sta- 
ble condition. We derived magnetic field strength, the ori- 
entation of the magnetic field, filling factor, and Dopplcr 
velocity using the Milnc-Eddington inversion code. In the 
following sections, we describe the details of the observa- 
tion in §2, analyze Stokes V area asymmetry in §3, show 
the statistical results of magnetic and Doppler fields of 
UDs in §4, and finally discuss and summarize our findings 
in §5. 

2. Observation 

The SP observation was performed from March 2 
through March 4, 2007, in parallel with the acquisition of 
the filtergram data analyzed in Kitai et al. (2007). The 
target was NOAA 10944 with a-typc sunspot in its de- 
caying phase. The sunspot invoked no flaring nor surging 
activity, and almost disintegrated on March 5. The region 
was located in the west side of the solar surface. The he- 
liocentric coordinate of NOAA 10944 was (S6°, W17°) on 
March 2, (S5°, W30°) on March 3, and (S6°, W43°) on 
March 4. For more information, please refer to Kitai et 
al. (2007). 

With the SP, a Normal Map mode was carried out from 
00:10 UT to 00:50 UT on the three consecutive days. The 
Normal Map mode scans an area with an integration time 
of 4.8 s per slit position. The observation covers the field- 
of-view (FOV) of 80"x80" with a polarimetric accuracy 
of 0.1%. The spatial pixel size was 0"159 in slit direction 
and 0"147 in step direction. The spectral FOV covers 
two absorption lines of Fe I 6301. 5A (g e ff =1.66) and Fe 
I 6302. 5A (g =2.5). The SSW routine spjprep.pro was ap- 
plied for the purpose of dark subtraction and flat fielding. 

3. Stokes V area asymmetry 

Figure 1 (top) shows the map of Stokes Varea asymme- 
try in Fe I 6302. 5 A line on March 2. The area asymmetry 
of Stokes V profiles provides an indication of a large gra- 
dient of field strength or line-of-sight velocity (Solanki & 
Stenflo 1984, Stenflo & Harvey 1985; Grossmann-Doerth 
et al. 1988, 1989; Sanchez Almeida & Lites 1992). We 
took the definition of the Stokes V area asymmetry (S A) 
as 
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I bl ue\y^-L d \V\d^ 

f b lue\V\d\ + f red \V\d\ 



(1) 



where each integration is performed over the area of the 
blue or red lobe of V profiles. In the dark core region 
(x=250", y=0"), SA is noisy, probably due to the blending 
of molecular lines at its cool temperature. Except for the 
dark core, the umbra has a negligible value of SA. This 
is because of the suppression of convective motions in the 
presence of strong magnetic field, which produces a small 
line-of-sight velocity in the umbra. Morinaga et al. (2007) 
reported a smaller asymmetry in the center of the pore 
than its surrounding, which is consistent with our result. 
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Fig. 1. Stokes V area asymmetry map (top) and contin- 
uum map at 6303A (bottom) on March 2. The white curves 
show the smoothed contours of the continuum intensity at 
Iquiet X 0.3. Here I qu i e t is the average intensity of the quiet 
region. The brightness of the continuum map is normalized 
by Iq U i e t- The Stokes profiles at the position (a) are shown 
in Fig. 2. 



4. Magnetic Structure and Doppler Velocity 
Distribution around UDs 

Inversion 

We applied a Milne-Eddington inversion code 
(Yokoyama et al. 2008, in preparation) to the Stokes 
spectra in. As the Stokes V profile in the umbra has a 
negligible area asymmetry (§3), it is reasonable to assume 
a Milne-Eddington atmosphere. The best-fit UD profiles 
at the position (a) indicated in figure 1 (bottom) are 
shown in figure 2. Figure 2 shows that Milnc-Eddington 
inversion produces excellent fits to the observations. 
The inversion code can derive 10 free parameters: three 
components of the magnetic field (strength, inclination, 
and azimuth), line-of-sight velocity, two parameters 
describing the linear dependence of the source function 
on optical depth, line strength, Dopplcr width, damping 
parameter, macroturbulcnt velocity, stray-light fraction, 
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Fig. 2. Observed (diamonds) and best-fit Stokes /, Q, U, 
V profiles (solid line) by the inversion at the position (a) 
indicated in Fig.l. The profiles are normalized by the Stokes 
/intensity at continuum wavelength (I c )- 



and a shift of the stray-light profile. The 180° ambiguity 
of the azimuth angle is determined by comparing with 
the potential fields calculated from the line-of-sight 
component of the magnetic field. The stray-light profile 
is the averaged Stokes / profile over the regions where the 
maximum polarization degree (p = sqrt(Q 2 + U 2 +V 2 )/I) 
along the line profile is larger than 0.2%. The stray-light 
represents the effect of a degradation of the polarization 
signals due to telescope diffraction and insufficient angu- 
lar resolution (Orozco Suarez et al. 2007). The magnetic 
filling factor, which represents the fraction of magnetized 
atmosphere, is computed as 1 — (stray- light fraction). 
The 2D maps and magnetic field vectors are converted 
to the local coordinate referring to the solar surface 
with the assumption that the solar surface is flat in our 
field of view, and, in the following, maps are presented 
as seen from the top. Dopplcr velocity is subtracted 
from the averaged Doppler velocity value inside the 
umbra. The observed Doppler velocity field includes 3 
min umbral oscillation, 5 min p-modc oscillation, and 
other instrumental effects. We did not filter out these 
effects, since our interest is the local variations around 
UDs, which is easily distinguishable from them (umbral 
oscillation, p-mode, etc.). 

We estimated the random error levels of the derived 
physical quantities from the standard deviation of the 
original map subtracted by a boxcar smoothed (width is 
0."3x0."3) map inside the umbra. The smoothed width 
(0."3x0."3) is chosen to be narrower than the typical UD 
size, to calculate the fluctuation level contribution from 
other sources. As a result, la error levels of field strength, 
field inclination, Doppler velocity, and filling factor are 13 
Gauss, 0.°7, 10 m s _1 , and 0.02, respectively. 

4-2. Identification of UDs 

To identify UDs, we took the image segmentation 
method explained in Sobotka et al. (1997). First, we made 
a boxcar smoothed (4x4 pixels) continuum map. Then, 
the original continuum map was divided by the smoothed 
one, that is, I c (original) / I c (smoothed). I c is the Stokes 
/ intensity at continuum wavelength. We set the empir- 
ical threshold value of 1.05 for identifying UDs, that is, 
an UD is more than 1.05 times brighter than its vicinity. 
To avoid the influence of the statistical noise, UDs whose 
areas are less than 3 pixels are excluded from the analysis. 
In this way, we identified 27 UDs on March 2 (Fig. 3), 35 
on March 3, and 25 on March 4. In total, 87 UDs are 
analyzed in this paper. 

In Kitai et al. (2007), UDs are classified into 3 cat- 
egories by their birth site, i.e., umbra, penumbra, and 
light-bridge origin. The three categories of UDs show dif- 
ferent characteristics with respect to their proper motion 
and temperature. However, we don't take care of the sub- 
classes in this paper, since we can not trace the birth site 
of each UD with our data sets. 
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Continuum intensity 



Fig. 3. Positions of 27 identified UDs on March 2. The di- 
rections of the arrows mean the horizontal orientation of mag- 
netic field at each position. The plots at the position (b) are 
shown in Fig. 4. 

4- 3. Results 

4- 3.1. Statistics over the three days 

We made spatial profiles of continuum intensity, field 
strength (B). field inclination (i), filling factor (J), and 
Doppler velocity (v) across each UD in the direction of 
the horizontal component of the magnetic field (shown in 
Fig. 3 with arrows). The plots for the UD (b) shown in 
figure 3 are shown in figure 4. The horizontal axis covers 
the spatial length of 2"2 (15 pixels), which is long enough 
to cover the entire UD. First, we decided the position of 
the UD (xjjd) and its background (xbq, 2 points at both 
side) by eye inspection of the local maximum and local 
minimums of the continuum intensity (Fig. 4 top). Second, 
we calculated AF = F(xjjd) — F(xbg), where F = B, i, 
/, or v. 

Figure 5 shows the histograms of the difference value 
(AF) of four physical parameters. In the statistical aver- 
age of all of the three day's data, UDs have relative blue 
shifts (Av average = -18ms _1 ). As for the magnetic field, 
however, the statistical averages (AB average = —7 Gauss, 



Ai 



average 



=0.°2, A/, 



average 



=0.007) are smaller than the 
error levels, i.e., 13 Gauss, 0.°7, and 0.02 respectively. 

In figure 6 we show four scatter plots (from top, field 
strength AB, field inclination Ai, filling factor Af, 
and Doppler velocity Av) against continuum intensity 
ratio UD/BG. Larger continuum intensity ratio means a 
brighter UD. Red (March 2), green (March 3), and blue 
(March 4) circles indicate the average values in intensity 
ratio bin A(UD/BG)=0.2 with error bars showing the 
standard deviation. There seems to be no correlation 
between continuum intensity ratio and Ai, and Af (the 
second and the third panels of Fig. 6). On the other 
hand, AB and Av have weak correlations with continuum 
intensity ratio (the top and the bottom panels of Fig. 6). 
That is, brighter UDs have weaker magnetic field, and 
larger blue shifts. 
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Fig. 4. Five plots around the UD indicated as (6) in figure 
3. From top to bottom, continuum intensity, field strength 
(£?), field inclination (i), filling factor (/), and Doppler veloc- 
ity (v) (positive means downflow). A thick bar below each 
plots means the span of pixels which satisfy UD's conditions 
stated in §4.2. xud = S,xgQ=5 and 10. -Byo = 2196Gauss, 
B B G = 2183Gauss, i UD =17.°0, i B G =16.°3, f UD = 0.880, 



= 0.888, v UD 



-70.1m s , vgQ = 35.6m s 
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Fig. 5. Histograms of UD—BG differences. The total sample 
number is 87. From top, field strength (AB), field inclination 
(Ai), filling factor (A/), and Doppler velocity (Ad). The 
black, gray, and light gray bars indicate the UDs on March 2, 
3, and 4, respectively. 
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Fig. 6. Scatter plots of field strength (AB), field inclination 
(Ai), filling factor (A/), and Doppler velocity (Aw) against 
continuum intensity ratio UD/BG. X signs are indicators of 
UDs on March 2 (red), March 3 (green), and March 4 (blue). 
The circles show the average values in intensity ratio bin 
A(UD/BG)=0.2 in each of the three observing days. The 
colored solid lines show the standard deviation error bars. 
The negative value of An means blue shift. 
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Table 1. Ccnter-to-limb variation of the averaged difference 
value UD-BG 



March 2 March 3 March 4 
S6°,W17° S5°,W30° S6°,W43° 




Field strength 


-17 


-6 
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(Gauss) AB 








Field inclination 


0.6 


0.1 


-0.1 


(degree) Ai 








Filling factor 


0.009 


0.004 


0.009 


A/ 








Doppler velocity 
(m s _1 ) Av 


-28 


-24 
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4-3.2. Center-to-limb variation 

As stated in §4.3.1, UD's magnetic field, compared 
to their surroundings, does not show distinct variations 
in the statistical average of all of the three days data. 
However, we found an interesting property in the daily 
statistics. The difference values AF show center-to-limb 
variation, as listed in table 1. |A£?|, |Ai|, and |Av| get 
smaller values as days go on, apart from A/. The sunspot 
was the closest to disk center on March 2. Observed near 
disk center, UDs show smaller field strength, larger field 
inclination (more horizontal), and relative blue shifts. On 
March 4, the sunspot was distant from the disk center. In 
this case, UDs and their BG show almost no difference in 
magnetic and Doppler field. As for the Doppler velocity, 
the projection effect may partly contribute to this result, 
since we detected the line-of-sight components. 

5. Discussion and Summary 

Our analysis revealed the magnetic structure and the 
velocity field distribution around UDs with seeing-frce, 
high-sensitive Hinode SP observation. The main results 
are as follows: 

1. The Stokes V profiles of UDs are virtually symmet- 
ric. 

2. In the statistical average of all of the three days 
data, UDs do not produce distinct variations from 
the surroundings in their magnetic conditions, while 
Doppler velocity shows effective blue shifts (—18 m 

3. The filling factor shows no difference on UDs. 

4. There are weak positive correlations between bright 
UDs and weaker field inclination, and between 
bright UDs and relative blue shifts. 

5. The local differences of field strength, field inclina- 
tion, and Doppler velocity on UDs show center-to- 
limb variation. 

5.1. Fe I formation height 

According to Tritschler & Schmidt (1997), the forma- 
tion height of the Fe I 6302. 5 A in umbrae is 180 km higher 
than the height of continuum optical depth re = 1 at 
500nm for the line core, and the formation height above 



Fig. 7. Schematic figure of the cusp-shaped UD magnetic 
field lines and the Fe I line formation height. The black solid 
lines are magnetic field lines. The central green part cor- 
responds to a weakly-magnetized hot gas, i.e., an UD. Blue 
and red arrows indicate upflow and downflow, respectively. 
The downflows, however, are an imaginary picture because 
we could not find them. The two gray thick lines indicate 
the line formation heights at disk center (lower) and at limb 
(upper). See text for further information. 

an UD is 130 km higher than tc — 1. In addition, the 
continuum level of the UD is found to be shifted to ^100 
km higher layers (Degenhardt & Lites 1993). But these 
values are estimated at disk center, where we can look 
into the deepest layer. When the target gets close to the 
limb, the formation height of spectral lines gets higher 
and higher. This may explain the center-to-limb varia- 
tion on UDs we've found in §4.3.2. The schematic view 
of an UD that accounts for the center-to-limb variation of 
magnetic and velocity fields is shown in figure 7. When 
we observe an UD at disk center, weakly-magnetized hot 
gas can reach the formation height of Fe I. The obser- 
vation reveals smaller field strength (Ai? = — 17Gauss), 
flatter field inclination (Ai =0.°6), and relative blue shifts 
(Av = —28m s" 1 ), though field inclination difference is less 
than the error fluctuation level. When we observe an UD 
far from disk center, the line formation height is higher 
than the UD occurring site. The observation supports 
this interpretation because we could not find large differ- 
ences between UDs and their BG on March 4. Our results 
directly support the model that the UDs are formed at 
deep photosphere, and have a cusp-shaped magnetic field. 
As for the downflows showed by red arrows in figure 7, 
the discussion is made in §5.3. 

There is a possibility that this center-to-limb variation 
is due to the evolutional phase difference of UDs in a 
decaying sunspot. No one has examined the change of 
UD characteristics in a developing, mature, and decaying 
sunspot as far as we know. Of course, the actual cause 
of this variation can be the mixture of the two, that is, 
the formation height difference and the evolutional phase 
difference. 
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5.2. filling factor 

As stated in §4.1, the filling factor is computed as 
1 — (stray-light fraction). In our inversion, the stray- light 
profile is the averaged Stokes / profile over the regions 
where the maximum polarization degree along the line 
profile is larger than 0.2%. According to the limb obser- 
vation of SP performed on March 16, 2007, the fraction of 
the scattered light in the continuum wavelength was 2% at 
a few arcsec away from the limb. On the other hand, the 
average stray-light fraction inside the umbra was ~10%. 
Thus the observed profile is considered to be composed 
of three components: unpolarized light due to telescope 
diffraction (<~2%), unpolarized light coming from the field- 
free atmosphere in the umbra, and polarized light coming 
from the magnetized gas. The filling factor indicates the 
fraction of the polarized light. 

Assume that UDs are the penetration of field-free hot 
gas (Socas-Navarro et al. 2004; Spruit & Scharmer 2006) 
and the penetration reaches the formation height of Fe 
I 6302A, the detection of small filling factor on UDs are 
predicted. However, we could not find a decrease of the 
filling factor on UDs (table 1 and Fig. 6). Higher resolu- 
tion observation is strongly needed for further discussion. 

5.3. Comparison with sunspot models 

5.3.1. The monolithic model 

A three dimensional simulation of UD phenomena was 
recently done by Schiissler & Voglcr (2006). They ex- 
plained UDs as a natural result of convection in a strong, 
initially monolithic magnetic field. In their simulation, 
most of the UDs had an elongated form with a central or 
threefold dark lane which separates the UD into two or 
three parts. At the end points of the dark lane there exist 
downflow patches. 

This picture is rather different from our observational 
impression. We observed almost circular UDs and no dark 
lanes inside them with Fe I continuum map. This is also 
the case with the blue/green continuum images obtained 
with Hinode SOT with spatial resolution 0."2 (Kitai et al. 
2007). However in recent days there has ever some evi- 
dence of dark lanes with high resolution data, for exam- 
ple, Bharti et al. (2007b) with G-band filtergrams using 
Hinode SOT. Thus far no one, including us, has succeeded 
in finding localized downflow patches at the end points of 
dark lanes (see Fig. 1 in Schiissler & Vogler 2006). There 
are two possibilities why we could not find localized down- 
flow patches: One is that, the downflow patches may be 
too small to be detected with Hinode's resolution limit. 
The other may be, because the continuum intensity (i.e., 
temperature) goes down when the gas flows downward, 
and the signal of the Stokes profile becomes too faint to 
be detected. 

We found another important phenomenon of UDs. 
It is already known empirically that some UDs occur 
and recur at the same location (Rimmclc 1997). With 
Hinode's secing-frce condition, we got the light curve of 
successive UDs for about 3 hours using green continuum 
data. Some UDs located in the center part of the umbra 




Fig. 8. Upper left: green continuum image on March 2, 2007. 
The position of a periodic UD is marked with a white square. 
Lower left: light curve of the center of the white square. Below 
the light curve shown diamond signs which indicate more than 
8% brighter than the surroundings. Right: Temporal change 
of zoom images of the white square region (pseudo-color dis- 
play). The time is written at upper left at each image in 
UT. 

showed oscillatory light curves. One example is shown 
in figure 8. The characteristic period of the oscillation 
is ~10 minutes, and the successive emergence of UDs 
continues over ^50 minutes. This intensity oscillation 
may be supporting evidence for the monolithic model, 
because this model indicates that oscillatory convection is 
the preferred mode in the first few Mm depth below the 
umbral photosphere by linear stability analysis (Weiss et 
al. 1990). 

5.3.2. The field- free intrusion model 

Another plausible model of UDs is the field-free intru- 
sion model. This model was proposed by Parker (1979), in 
which he discussed that the magnetic field lines beneath 
the umbra are divided into many separate flux tube bun- 
dles, like "spaghetti". The UDs would be a manifestation 
of field-free hot gas intrusion from below through the gaps 
of nonuniform magnetic field. 

A numerical study of the field-free intrusion model was 
done by Spruit & Scharmer (2006) for pcnumbral grains 
and an observational study was done by Rimmele (2008). 
They predict cusp-shaped magnetic field lines, smaller 
magnetic field, and upflow within the cusp. These charac- 
teristics of UDs are almost the same as those of the mono- 
lithic model, apart from the localized downflow patches. 
Field-free intrusion model also predicts downflow based 
on radiative cooling, but it's not concentrated on local- 
ized patches. The weak correlation between dark UDs and 
downward motion, shown in figure 6 may be evidence for 
downflow by radiative cooling. As for the oscillatory light 
curve (Fig. 8), however, the field-free intrusion model fails 
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to predict such oscillation. It would be very helpful if a fu- 
ture numerical simulation can make clear the possibility of 
brightness oscillation in the field-free hot gas surrounded 
by strong umbral magnetic field. 

As was discussed in this paper, the oscillatory bright- 
ening of UDs seems to be a key phenomenon to reach 
a conclusion on the origin of UDs. Detailed numerical 
studies and higher resolution observational, including 
the spectroscopic study of its temporal evolution on this 
phenomenon, are strongly needed in the near future. 

We are certain that we found many meaningful observa- 
tional properties of magnetic field around UDs. The rea- 
son why we succeed in deriving good correlations between 
UDs and the components of magnetic field, owes greatly 
to the really stable, sensitive, and high-resolutional obser- 
vation performed by spcctro-polarimetry on board of the 
Hinode Solar Optical Telescope. Hinodc is a Japanese mis- 
sion developed and launched by ISAS/JAXA, with NAOJ 
as domestic partner and NASA and STFC (UK) as in- 
ternational partners. It is operated by these agencies in 
co-operation with ESA and NSC (Norway). 
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Sports, Science and Technology (MEXT) of Japan, and 
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